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Abstract 

By performing density functional theory (DFT) +U calculations, we systematically study the 
electronic, mechanical, tensile, and thermodynamic properties of Am02- The experimentally ob- 
served antiferromagnetic insulating feature [J. Chem. Phys. 63, 3174 (1975)] is successfully repro- 
duced. It is found that the chemical bonding character in Am02 is similar to that in PuC>2, with 
smaller charge transfer and stronger covalent interactions between americium and oxygen atoms. 
The valence band maximum and conduction band minimum are contributed by 2p-5f hybridized 
and 5/ electronic states respectively. The elastic constants and various moduli are calculated, 
which show that A111O2 is less stable against shear forces than PuC>2. The stress-strain relation- 
ship of Am02 is examined along the three low-index directions by employing the first-principles 
computational tensile test method. It is found that similar to Pu02, the [100] and [111] directions 
are the strongest and weakest tensile directions, respectively, but the theoretical tensile strengths 
of A111O2 are smaller than those of Pu02- The phonon dispersion curves of A111O2 are calculated 
and the heat capacities as well as lattice expansion curve are subsequently determined. The lattice 
thermal conductance of Am02 is further evaluated and compared with attainable experiments. 
Our present work integrally reveals various physical properties of A111O2 and can be referenced for 
technological applications of A111O2 based materials. 
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I. INTRODUCTION 



Actinide based materials possess interesting physical behaviors due to the existence of 5/ 
electrons and have attracted extensive attentions jl-7], among which the actinide dioxides 
(An0 2 ) have been investigated most widely. At their ground state, all actinide dioxides 
crystallize in face-centered-cubic (CaF2-like) structure. In this arrangement, each actinide 
atom, located at the center of an oxygen cube, is expected in the ionic limit to yield four of 
its electrons to the surrounding oxygen atoms. This would lead to formal integer 5/ orbital 
populations ranging from f° (Th) to f 10 (Fm). Among the dioxide series, Pu0 2 and Am02 
occupy the intermediate zone, before Pu0 2 (Th0 2 — > Np0 2 ) the 5/ electronic states appear 
as localized states above the oxygen-2p states, while after Am0 2 (Cm0 2 — > Fm0 2 ) the 5/ 
electronic states appear as localized states below the oxygen- 2p states. Thus the actinide 



5/ electronic states in Am0 2 as well as in Pu0 2 have the largest over 
2p states. However, different from the widespread researches on Pu0 2 



ap with oxygen- 



7K12J, Am0 2 have 



received very few theoretical concerns until now. Therefore following our previous systematic 



111 ], here 



investigations on the ground-state properties and high-pressure behaviors of Pu0 2 
we further investigate the corresponding physical properties of Am0 2 . 

Another reason for a thorough study of Am0 2 is its important role played in the latest 
nuclear reactor fuels. During the burning cycle of U0 2 in conventional fission nuclear reac- 
tors, some considerable amounts of plutonium and neptunium isotopes, as well as smaller 
quantities of minor actinides such as Am, Cm, Bk, and Cf, emerge in the reaction waste 13]. 
Since such radioactive waste results in very troublesome long-term storage requirements, it is 
then suggested that the produced transuranium element dioxides could be reprocessed from 



the burnt fuel and used as alternative fuel in a new generation of reactors 



Ij4l7j. Therefore, 



the mechanical and thermodynamic properties of these dioxides also need to be investigated 



ll| and Np0 2 (l8|, 



and presented. Similar to our previous theoretical studies on Pu0 2 
here we will systematically obtain various physical properties of Am0 2 and compare them 
with attainable experiments as well as possible. 

Different physical properties of Am0 2 have been investigated in experiments ever since 
1969 



19 



20J . The earlier magnetic susceptibility measurements reveal that Am0 2 undergoes 



antiferromagnetic transition at 8.5±0.5 K 



20j. Recently, along with the development of 



computing abilities and numerical methods for applying quantum mechanics calculations, 
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more and more theoretical researches are now being carried out on multi-phase properties of 
different actinide dioxides. Although problems still exist like on the ground-state magnetic 



orders 



21 



251 ] . many theoretically obtained physical properties of dioxides like U0 2 and 



PuC>2 already accord very well with experiments. For example, the optical spectrum of UO2 
given by employing the density functional theory (DFT) +U method 26] is well consistent 
with the optical reflectance measurement 27], which has inspired experimentalist 28] to 



choose the DFT+[7 results of PuC>2 



26] as a criteria in determining the melting temperature 



of this material. Also remarkably, the DFT+[7 predicted phonon spectrum of PuC>2 11] 
has been confirmed by the most recent experimental report [29|[. Based on these previous 
successes, in the present paper we intend to reveal the ground-state properties for AH1O2, 
including electronic, mechanical, and thermodynamic aspects, as well as its high-pressure 
behaviors. 



II. COMPUTATION METHODS 



Our DFT calculations are carried out using the Vienna ab initio simulations package 



31] . The ex- 



(VASP) [30] with the projected-augmented-wave (PAW) potential methods 
change and correlation effects are described by local density approximation (LDA) and gen- 



eralized gradient approximation (GGA) in the Perdew-Burke-Ernzerhof (PBE) form [32 ]. 
The plane-wave basis set is limited by an energy cutoff of 500 eV. Integrations over the 
Brillouin Zone of the cubic cell containing 4 americium and 8 oxygen atoms are done on 



a 9x9x9 k points mesh generated by the Monkhorst-Pack [33j method, which is sufficient 
for an energy convergence of less than l.OxlO -4 eV per atom. The strong on-site Coulomb 
repulsion among the localized americium 5/ electrons is described by the DFT+U formal- 



ism formulated by Dudarev et al. 34J-I36] . In this method, only the difference between the 
spherically averaged screened Coulomb energy U and the exchange energy J is significant 
for the total LDA/GGA energy functional. And in this paper, the Coulomb U is treated 
as a variable, while the exchange energy is set to be a constant J=0.75 eV. Since only the 
difference between U and J is significant, thus, we will henceforth label them as one single 
parameter, labeled as U for simplicity, while keeping in mind that the nonzero J has been 
used during calculations. 

In order to calculate the elastic constants for A111O2, we enforced small strains on the 
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equilibrium cubic cell. For small strain e, Hooke's law is valid and the crystal energy E(V, e) 



can be expanded SIS db Taylor series 



6 6 



E(V, e) = E(V , 0) + V J2^ + yJ2J2 C ^ e i + °({ e * »' C 1 ) 

i=l i=l j=X 

where E(Vq, 0) is the total energy at the equilibrium volume Vo of the cell without strains, 
Cij are the elastic constants, and e is the strain tensor which has matrix elements (i,j = l, 
2, and 3) defined as 



"13 



\ 



ei §e 6 |e 5 ~e 6 e 2 |e 4 |e 5 



e 2 |e 4 ~e 5 |e 4 e 3 



|e 5 |e 4 e 3 



(2) 



V / 

Note that we have used the Voigt notation in the equation above which replaces xx, yy, zz, 
yz, xz, and xy by 1, 2, 3, 4, 5, and 6, respectively. For cubic structure, there are only three 
independent elastic constants, Cn, Cu, and C 44 . So, we can employ three different strains 
to calculate them as follows, e 1 = (<5,5,0,0,0,0), e 2 =(5,-<5,0,0,0,0), e 3 =(0,0,0,<5,<5,0). The strain 
amplitude 5 is varied in steps of 0.02 from —0.06 to 0.06 and the total energies E(V, S) at 
these strain steps are calculated. 

For semiconductors, the Helmholtz free energy F can be expressed as 

F(V,T) = E(V) + F mb (V,T), (3) 

where E(V) stands for the ground-state electronic energy, and F v n,(V,T) is the phonon free 
energy at a given unit cell volume V. Under quasi-hamonic approximation (QHA), F vi fy(V, T) 
can be evaluated by 

F vib (V,T) = k B TY,^ 
i>q 

where w 3 -(q, V) is the phonon frequency of the jth phonon mode with wave vector q at fixed 
V, and ks is the Boltzmann constant. The total specific heat of the crystal is the sum of all 



2smh(W)) 
V 2k B T J 
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phonon modes over the Brillouin zone (BZ), 

C B (T) = X)c B j(q,T). (5) 

3,1 

c v j (q, T) is the mode contribution to the specific heat defined as 

c „ j(q , t) = £ (5^) sinh2[ ^ (q V )/2fer] - <«) 



The mode Griineisen parameter 7j(q) describing the phonon frequency shift with respect to 
the volume can be calculated by 

_ #n^(q,U)] 

The lattice thermal conductivity k for a material can be written as a sum over one 



longitudinal (kla) and two transverse (kta and kta') acoustic phonon branches }38l l3Qf| . 



K= K LA + K TA + K TA ,. (8) 

These partial thermal conductivities can be calculated differently, depending on the specific 
mechanisms for phonon scattering rates l/r c , where r c is the relaxation time. At relative 
high temperatures, the dominant mechanism for phonon scattering is the normal and Umk- 
lapp phonon-phonon processes (l/r c = 1/tjv + m which mainly the acoustic phonon 

branches interact with each other and transport heat. Using the Debye-Callaway model 
3^, the partial conductivities k { (i corresponds to TA, TA', or LA modes) can be 
expressed as 

1 2 



rQi/T rl{x)x^ , 
Jo T i N {e*-l) 2U '- L 



(9) 



3 Wo (e x - l) 2 r 0l/T Z&3*fg w r 

I Jo t^t^ 

where Oj is the longitudinal (transverse) Debye temperature, x = hw/ksT, and Ci = 
k B /2ir 2 h 3 Vi. Here, h is the Plank constant, k B is the Boltzmann constant, u is the phonon 
frequency, and V{ is the longitudinal or transverse acoustic phonon velocity. The normal 
phonon scattering and Umklapp phonon-phonon scattering are written as 

i kW LA v fksV^ {10) 



r^ A (x) Mh 2 v\ A V h 

1 _ ^bItA/TA'^ k B 



r T N AITA \x) Mh^ TA/TA , h 



xT\ (11) 
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TABLE I: Calculated lattice parameters (ao), bulk modulus (Bo), pressure derivative of bulk 
modulus (B' ), spin moment of each americium ion (/x mag ), and energy band gap (E g ) for antifer- 
romagnetic A111O2 by different calculational methods. Other experimental and theoretical results 
are also included for comparisons. 



Method 


ao 


B 


B o 


f^mag 


Eg 




(A) 


(GPa) 




(Pb) 


(eV) 


LDA+C7 ([7=0) 


5.280 


217.7 


4.53 


4.56 





LDA+C7 ([7=4) 


5.351 


189.6 


4.71 


4.94 


0.7 


GGA+[7 (17=0) 


5.401 


179.3 


4.10 


4.77 





GGA+C7 (U=4) 


5.484 


140.1 


4.81 


5.15 


1.0 


Expt. 


5.383 a ,5.375 fe 










HSE 


5.37 c 






5.1 c 


1.6 C 



a Reference 42 . 
b Reference 43. 
c Reference |7. 



and 

1 h-f 2 (k B x 2 



^ x 2 T 3 e" e ^ 3T , 



(12) 



respectively, where 7 is the Griineisen parameter, V is the volume per atom, and M is 
the average mass of an atom in the crystal. With reasonable expressions of the Debye 
temperature and acoustic Griineisen parameter to describe the harmonic phonon branches 
and the anharmonic interactions between different phonon branches, Eq. (Q can provide 
reasonable predictions for a material's thermal conductivity. 



III. RESULTS AND DISCUSSIONS 

A. Structural and electronic properties 

Firstly, the equilibrium lattice parameter of AH1O2 is determined from total-energy cal- 
culations by the DFT+U methods. The total energies are calculated in a wide range of 



6 



lattice constants and fitted to the Brich-Murnaghan equation of states (EOS) 41] 

-l 3 



E(V) = E + 



9V B 



+ 



16 

V 



V \ 2/3 

— I - 1 

V 



B' 



2 r 



(13) 



where E is the electronic energy of Am02 at its equilibrium volume Vq, B is the bulk 
modulus, and B' is the pressure derivative of the bulk modulus. The numerically fitted 
lattice constants, B Q , and B' are summarized in Table I. For comparison, the experimental 



values 



42 



43| and HSE theoretical results jz| are also listed. One can see that the LDA and 



LDA+[7 results on the lattice constant a® are both smaller, while the GGA and GGA+U 
results are both larger than the experimental values. This fact reflects the overbinding effect 
of LDA and the underbinding effect of GGA on the chemical bonding strengths between 
different atoms. However, all the obtained lattice constants are not too far away from the 
experimental values. As for bulk modulus Bo, it displays strong dependence on the U values 
for both LDA+[7 and GGA+U approaches. As U increases from to 4 eV, B decreases 
monotonously by 12.9% and 21.9% for LDA+?7 and GGA+U calculations, respectively. Due 
to overbinding effect of LDA and underbinding effect of GGA, the LDA+U results of B are 
always higher than the GGA+U ones. We have also calculated the bulk modulus Bo from 
empirical formula based on the elastic constants (as shown in Table III), which turns out to 
be very close to the result by EOS fitting, indicating that our calculations are self-consistent. 
Moreover, through EOS fitting, we also systematically obtain the pressure derivative of bulk 
modulus by using different methods. The corresponding results are also listed in Table I. 

Through electronic structure calculations, we find that without employing the +U method 
describing the strong electron correlation effect, A111O2 is wrongly described to be metallic 
in its ground state. Figure 1 shows the projected density of states (PDOS) for A111O2 
obtained with different computational methods. Since spin-orbital coupling (SOC) effects 
might be important to electronic structure descriptions of actinide dioxides, here we consider 
SOC in all calculations. In the LDA/GGA calculational level, SOC makes americium 5/ 
orbitals to split into the lower j—5/2 (^5/2) and higher j—7/2 (^7/2) orbitals, occupying below 
and above the Fermi energy respectively, as shown in Figs. 1(a) and 1(c). The obtained 
electronic states show metallic characters with considerable / 5 / 2 states distributing at the 
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FIG. 1: (Color online). Electronic density of states for A111O2 in its AFM state calculated by 
using the (a) LDA, (b) LDA+U, (c) GGA, and (d) GGA+U methods. The spin-orbit coupling is 
included in all calculations. The value for the U parameter is chosen to be 4 eV for all LDA+C/ 
and GGA+U calculations. The Fermi energy is set at zero. 

Fermi energy. However, after turning on the [/-parameter to describe the strong on-site 
electron correlation effects of americium 5/ electrons, the state is further split into two 
separate sub-bands, distributing at the two sides of the Fermi energy. As shown in Figs. 
1(b) and 1(d), an obvious energy band gap appears in LDA/GGA+f/ calculations. Within 
LDA and GGA+U calculations, the energy band gap of Am02 is found to be 0.7 and 1.0 
eV respectively, as shown in Table I. As will be discussed below, the difference between the 
obtained energy band gaps might come from the different hybridization strengths of the 
americium 5/ and 6d states in LDA/GGA+?/ calculations. The obtained insulating ground 
state of AH1O2 is in agreement with previous hybrid density functional result |7|], where the 
energy band gap is found to be 1.6 eV. From the obtained PDOS distributions shown in 
Figs. 1(b) and 1(d), we can see that the americium 5/ and 6d states hybridize below the 
Fermi energies, and this hybridization is stronger in GGA+U than in LDA+U calculations. 
The stronger hybridization moves the 5/ states lower, and results in the fact that the band 
gap is larger in GGA+U calculations. The calculated local spin moment of each americium 
ion are 4.94 fis and 5.15 \xb in the LDA+£7 and GGA+U calculations respectively, which 
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FIG. 2: (Color online). The cubic cell of A111O2 with Fm3m space group, and a contour map 
depicting the charge density distribution inside the (110) plane. Large cyan and small white 
balls represent for americium and oxygen atoms respectively. The charge density distribution is 
calculated by using the GGA+C/ method with the U parameter of 4 eV. 

are both similar to the 5.1 \lb result obtained in hybrid density functional calculations 

To further understand the chemical bonding characters of AmC>2, we present in Fig. 
2 the charge density distribution inside the (110) plane. Similar to what we have found 



in our previous studies on other actinide dioxides [111, [18], the charge density distribution 
around the Am and O ions are both spherical, with small deformations in the directions 
toward their nearest neighboring ions. This character is typical for mixed ionic /covalent 



chemical bondings in actinide dioxides By employing the Bader analysis [44j, |45| which 
decomposes the charge density distribution in real space following topological methods, we 
further calculate the atomic charges for Am and O ions and present the results in Table II. 
Note that although we have included the core charge in charge density calculations, only 
the valence charge results are listed. Given that the valence electron numbers of Am and 
O atoms are 17 and 6, about 2.26 electrons transfer from each Am to O atoms. In our 
previous studies, we have also applied similar DFT+U calculations for UO2, NpC>2, and 
PuO"2 [ill, [3] and found that the charge transfer from each actinide to oxygen atoms is 2.56 
for UO2, 2.48 for NpO"2, and 2.40 for PuO^. From the charge transfer results, we can see that 
with increasing the atomic number of the actinide elements, the chemical bondings between 
actinide and oxygen atoms show less ionic characters. In addition, the bonding lengths of 
the dioxides are also listed in Table II. Accompanying with the weakening of ionicity, the 
bonding lengths between actinide and oxygen atoms decrease with increasing the atomic 



TABLE II: Calculated atomic charges according to Bader partitioning as well as the bonding 



lengths of different actinide dioxides. 





Qb(A) 


Qb(O) 


bonding length 


Compound 


(e) 


(e) 


(A) 


Am0 2 


14.74 


7.13 


2.27 


Pu0 2 a 


13.60 


7.20 


2.32 


Np0 2 fe 


12.52 


7.24 


2.34 


U0 2 a 


11.44 


7.28 


2.36 


Th0 2 c 


9.34 


7.33 


2.43 



a Reference 11 . 
b Reference 18. 
c Reference 146. 



number of actinide elements. 



B. Mechanical properties 

Based on the Hooke's law, we systematically calculate the three independent elastic 
constants for Am0 2 within the GGA+U formalism. The results are listed in Table III 
together with our previous results on other actinide dioxide materials. At the value of U=4 
eV, the calculated values of Cn, Ci 2 , and C44 are 250.4, 87.0, and 55.3 GPa respectively. 
Mechanically, the phase of Am0 2 is stable due to the fact that its elastic constants satisfy 



the following mechanical stability criteria [37| of cubic structure: 

C n > 0, C 44 > 0, C n > |C 12 |, (C n + 2C 12 ) > 0. (14) 

Among the actinide dioxide series listed in Table III, the elastic constants of Am0 2 have the 
lowest values. Based on the elastic constants, we can further calculate the bulk and shear 
moduli from the Voigt-Reuss-Hill (VRH) approximations. Firstly, the Voigt and Reuss 
limits on the bulk (By and Br) and shear moduli (Gy and Gr) are calculated following 



the expressions in Refs. |47j, |48|, then the bulk and shear moduli can be approximately 
evaluated by B = \{By + Br) and G = \{Gy + Gr) 49]. In addition, the Young's modulus 

10 



E and Poisson's ratio v can also be evaluated from the elastic constants by E = %b+g 
and v = 2(3B+G) ' res P ec tively. Our results for mechanical properties of Am02 by using 
the GGA+U method are collected in Table III. The value of 141.5 GPa for bulk modulus 
B is in agreement with the EOS-fitting result of 140.1 GPa. Besides, in comparison with 
other actinide dioxides included in Table III, the value of B for Am02 is much lower. This 
means that when external pressures are applied on these actinide dioxides, AH1O2 will be 
the easiest to get cracked. For the shear modulus, we can see that PuO"2 and AmO"2 have 
relatively smaller values than UO2, NpO"2, and Th02- This result comes from the fact that 
the oxygen 2p and actinide 5/ electronic states have more overlaps in PuC>2 and AH1O2, and 
subsequently the chemical bondings between oxygen and actinide atoms show more covalent 
characters in PuO"2 and Am02- Considering that the values of B and G represent for the 
ability to resist external pressures and shearing forces respectively, Pugh et al. su gge sted 
that their quotient B/G can be used to scale the ductility or brittleness of a solid |50|] . A 
high value of B/G means that the material has more ductility than brittleness, and a critic 
value of 1.75 is suggested to distinguish ductile or brittle materials. Our calculated value 
of 2.19 for B/G means that Am02 is a ductile material. At last, our calculated Young's 
modulus and Poisson's ratio for AH1O2 are 168.4 GPa and 0.302 respectively, both of which 
are among the typical values of insulating oxide materials. 

The Debye temperature (O) of a solid can be determined following the Debye theory, 
which describes atomic vibrations to be elastic waves, and subsequently O is related to the 



sound velocity of a material |5l|. Specifically 



o h ( 3n V' 3 

% = T B \^a) "™ (15) 

where h and fee are Planck and Boltzmann constants, respectively, n is the number of 
atoms in one unit cell, Q is the volume per unit cell, and v m is the averaged sound velocity. 
Approximately, the averaged sound velocity can be calculated by 



1/2 1 
3' + 



1-1/3 

(16) 



where v t = ■sjG/ p (p is the density of a material) and vi = \/ (3B + AG)/3p are the transverse 
and longitudinal elastic wave velocities, B and G are the bulk and shear moduli. Our 
calculated sound velocities and Debye temperature of A111O2 are also listed in Table III. 
For comparison, our previous results for other actinide dioxides are also included. We 
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TABLE III: Calculated elastic constants, various moduli, and Poisson's ratio v for Am0 2 . For 
comparisons, previous experimental values and other theoretical results are also listed. 







Cl2 


C44 


B 


G 


E 


V 




vi 


Vm 


e 




(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 


(GPa) 




(m/s) 


(m/s) 


(m/s) 


(K) 


Am0 2 


250.4 


87.0 


55.3 


141.5 


64.7 


168.4 


0.302 


2416.6 


4534.1 


2699.8 


335.7 


Pu0 2 a 


256.5 


167.9 


59.2 


197 


52.7 


145.2 


0.377 


2477.6 


5206.3 


2787.0 


354.5 


Np0 2 6 


363.6 


118.8 


57.4 


200 


78.1 


207.5 


0.327 


2835.0 


5566.5 


3176.8 


401.2 


uo 2 a 


389.3 


138.9 


71.3 


222.4 


89.5 


236.7 


0.323 


2841.8 


5552.7 


3183.4 


398.1 


Th0 2 c 


349.5 


111.4 


70.6 


191 


87.1 


226.8 


0.302 


2969.1 


5575.5 


3317.3 


402.6 



a Reference 11 . 
b Reference 18. 
c Reference 46. 



can see from Table III that the v t and vi values of Am0 2 are lower than those of Pu0 2 , 
U0 2 , Np0 2 , and Th0 2 . Since the sound velocity of a compound is always related to its 
lattice thermal conductivity, this result reflects that Am0 2 has smaller thermal transport 
properties. Because the Debye temperature is proportional to v m , Am0 2 among the listed 
actinide dioxides possesses the lowest Debye temperature. 



C. Theoretical tensile strength 

Based on the obtained atomic structure of Am0 2 , we further study its theoretical tensile 
strength. The ideal strength of materials is the stress that is required to force deformation 
or fracture at the elastic instability. Although the strength of a real crystal can be changed 
by the existing cracks, dislocations, grain boundaries, and other microstructural features, its 
theoretical value can never be raised, i.e., the theoretical strength sets an upper bound on the 



attainable stress. Here, we employ a first-principles computational tensile test (FPCTT) 52] 
to calculate the stress-strain relationship and obtain the ideal tensile strength by deforming 
the Am0 2 crystals to failure. The anisotropy of the tensile strength is tested by pulling the 
initial fluorite structure along the low-index [001], [110], and [111] directions. As shown in 
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FIG. 3: (Color online). Schematic illustration of structure models for applying tensions along the 
[100] (a), [110] (b), and [111] (c) directions, (d) and (e) Total energies per formula unit and stresses 
for AH1O2 being pulled along the [100], [110], and [111] directions. 

Figs. 3(a)-3(c), three geometric structures are constructed to investigate the tensile strengths 
in the three high-symmetry crystallographic directions: Fig. 3(a) shows a general fluorite 
structure of Am02 with four Am and eight O atoms; Fig. 3(b) a body-centered tetragonal 
unit cell with two Am and four O atoms; and Fig. 3(c) an orthorhombic unit cell with six 
Am and twelve O atoms. The tensile stress is calculated according to the Nielsen-Martin 
sdleme = where e„ is the strain t en S0 , («, ,=1, 2, 3, m d „ is t ee votane 

at the given tensile strain. Tensile process along the [001], [110], and [111] directions is 
implemented by increasing the lattice constants along these three orientations step by step. 
At each step, the structure is fully relaxed until all other five stress components vanished 
except that in the tensile direction. 

The calculated total energies and stresses as functions of uniaxial tensile strains for AH1O2 
along the [100], [110], and [111] directions are shown in Figs. 3(d) and 3(e). One can see 
that with increasing the tensile strains along all the three directions, the energies increase 
with inflexions. By differentiating the energy-strain curves, we find the inflexion points 
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TABLE IV: Calculated strain energy (Estrain} and stress maxima at the critic points in the tensile 



processes. 





Strain 


Egtrain 


Stress 


Direction 




(eV/atom) 


(GPa) 


[100] 


0.42 


2.87 


51.8 


[110] 


0.24 


0.64 


13.4 


[111] 


0.18 


0.35 


4.6 



along the [100], [110], and [111] directions to be the strain values of 0.42, 0.24, and 0.18 
respectively. At these strain values, the corresponding energy increases are 2.87, 0.64, and 
0.35 eV per atom in A111O2, and the corresponding stress maxima are 51.8, 13.4, and 7.6 GPa 
respectively, as shown in Table IV. Since the stress-strain curves change from upward into 
downward at these inflexion points, the stress maxima values are regarded as the theoretical 
tensile strengths. Clearly along the three low-index directions, the [100] direction has the 
largest, while the [111] direction has the smallest tensile strength. This result is similar to 



what we have found in another actinide dioxide PuO"2 [11] . We can also understand this 
result from chemical bonding analysis. Along the [100] direction, there are eight Am-0 
bonds per formula unit for fluorite AmO"2, and the angle of all eight bonds with respect to 
the pulling direction is 45°. In comparison, along the [110] and [111] directions only four 
and two Am-0 bonds make an angle of 45° with the pulling direction. The other four bonds 
along the [110] direction are vertical to the pulling direction, while the other six bonds make 
an angle of about 70.5° with the pulling direction. It is evident that the bonds vertical to 
the pulling direction have no contributions to the tensile strength, and the bonds parallel to 
the pulling direction are easy to fracture under tensile deformation. Therefore, the tensile 
strength along the [100] direction is stronger than that along other two directions. Along 
all the three direction, no structural transition has been observed under external strains in 
our present study. Besides, we note that the stress in [110] direction experiences an abrupt 
decrease process after strain up to 0.24. This is due to the fact that the corresponding four 
Am-0 bonds making an angle of 45° with the pulling direction have been pulled to fracture 
under that strain. In comparison with PuC>2 11], we can see that the theoretical tensile 
strengths of AmO"2 are much smaller along all the three low-index directions, indicating that 
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Am02 is easier to get deformed during tensile processes. 

D. Phonon dispersion curves and thermodynamic properties 

To our knowledge, no experimental or theoretical phonon frequency results have been 
published for Am02- Here we carry out the first theoretical try in obtaining its phonon 
dispersion curves, as well as thermodynamic properties. The density functional perturbation 
theory (DFPT) and a 2x2x2 supercell containing 96 atoms is adopted to calculate the 
force constants of Am0 2 . Since our previously obtained phonon dispersion curves result by 
using the same methods for PuC>2 accords moderately well with subsequent experimental 
measurements, we are confident that our present results for A111O2 can be effective to certain 
degrees. The Born effective charges are firstly calculated to reveal the dielectric properties 
of AH1O2, and to modify the LO-TO (longitudinal optical and transverse optical branches 
of phonon) splitting effects. For fluorite AE1O2, the effective charge tensors for both Am 
and O are isotropic because of their position symmetry. After calculation, the Born effective 
charges of Am and O ions are found to be Z* ArD = 4.51 and Z* = —2.26, respectively, within 
GGA+U formalism with the choice of £7=4.0 eV. We present in the left panel of Fig. 4 
(solid lines) our calculated phonon dispersion curves for AmCV In the fluorite primitive 
cell, there are three atoms (one Am and two O atoms). Therefore, nine branches of phonon 
dispersion curves exist. One can see that the LO-TO splitting at T point becomes evident 
by the inclusion of polarization effects. In addition, due to the fact that americium atom is 
heavier than oxygen atom, the vibration frequency of americium atom is lower than that of 
oxygen atom. As a consequence, the phonon density of states for Am02 can be viewed as 
two parts, as shown in the right panel of Fig. 4. One is the part lower than 4.0 THz where 
the main contribution comes from the americium sublattice while the other part higher than 
4.0 THz are dominated by the dynamics of the light oxygen atoms. 

To investigate the effects of strains on the phonon dispersion, we also calculate the phonon 
spectrum for Am0 2 that is stretched with an isotropic strain of 0.03, and show the results 
in Fig. 4 with dashed lines. Clearly after being expanded, both the acoustic and optical 
branches of phonon experience frequency shift downward, indicating the weakening of atomic 
interactions. 

The thermodynamic properties of a material are connected to its phonon dispersion 
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k r x w r l pDOS 



FIG. 4: (Color online). Phonon dispersion curves along the high symmetry lines of the Brillouin 
Zone (left panel) and phonon density of states (right panel) for A111O2 at the equilibrium and 
expanded volumes. Solid and dashed lines represent for the results at the equilibrium and expanded 
volumes respectively. The expansion ratio is 0.03 for the expanded lattice calculations. 
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FIG. 5: (Color online). Temperature dependence of our calculated linear thermal expansion (solid 
line) and bulk modulus ratio (dashed line), and the experimentally measured lattice expansion 
results for A111O2. 



curves. By systematically calculating the Helmholtz free energies at different lattice con- 
stants and different temperatures, we can determine the lowest-energy lattice constant for 
AH1O2 at different temperatures. As expected, the lattice constant enlarges as temperature 
increases. And the obtained lattice expansion curve for Am02 is shown in Fig. 5, together 
with the experimental results by Taylar et al. 54j and Fahey et al. 55]. Similar linear 
relationships can be found in our theoretical and previous experimental results. The small 
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FIG. 6: (Color online). Calculated heat capacities of Am02 at constant volume (C v ) and constant 
pressure (C p ), and corresponding experimental results as well as other theoretical results. 

difference on the increasing ratios comes from the fact that we only consider harmonic vi- 
brations during calculations of the phonon dispersion curves and the Helmholtz free energy. 
By fitting the Brich-Murnaghan EOS for Am02 at different temperatures for Am02, we 
also obtain the temperature dependence of its bulk modulus B. And our result of B/Bq as 
a function of temperature T is also shown in Fig. 5. One can see that as temperature in- 
creases, the bulk modulus decreases, indicating that Am02 is easier to get cracked at higher 
temperatures. 

Within QHA the heat capacity at constant volumes can be calculated following Eqs. (J5j) 
and Besides, the heat capacity at constant pressure C p can be further evaluated by 
using the relationship 

C p -C v = a 2 v (T)B(T)V(T)T. 

We then calculate the heat capacities for AmO ? , and show the results in Fig. 6. For compar- 
ison, the experimental data by Thiriet et al. 56] and Nishi et al. [43] are also shown in Fig. 
6. We can see that with increasing the temperature, the value of C p increases continuously, 
while the value of C v approaches to a constant of 3R (R is the gas constant). As clearly 
shown, our theoretical result for C p is in good agreement with the experimental measure- 



ments by Thiriet et al. 



56| at the temperature range from 300 to 500 K. As temperature 



further increases, the discrepancy between our result and the experimental values 43|, [56] 
become larger. We conclude that the disagreement of C p between theory and experiments in 
the high temperature range mainly originates from the QHA we used. Due to anharmonic 
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effects, our obtained lattice expansion curve is lower than experimental results, as shown in 
Fig. 5. The underestimation on the lattice expansion of Am02 finally results in the discrep- 
ancies between our calculated C p values and corresponding experimental results through the 
above relationship between C p and C v . 



E. Griineisen parameters and thermal conductivity 

Based on the calculated phonon dispersion curves of Am02 at both the equilibrium 
and expanded lattice volumes, we can further calculate the Griineisen parameters for each 
vibration mode of Am02- Our obtained Griineisen parameters at different k points are shown 
in Fig. 7. Due to the fact that after lattice expansion, atomic interactions are weakened and 
phonon frequencies are lowered, the Griineisen parameters are positive for almost all phonon 
branches. The only exception lies around the W point, where the Griineisen parameter for 
the second acoustic branch of phonon has a small negative value. This result comes from 
the coupling behavior of the two lowest-frequency acoustic phonon branches at the W point. 
After lattice expansion, the coupling disappears and the dispersion curves of these two 
branches separate with each other, causing that the frequency of the higher branch moves 
upward. As a result, the Griineisen parameters for this phonon branch have small negative 
values around the W point. 

We have explained above that the Umklapp scattering ^ is proportional to 7 2 , thus, the 
Griineisen parameter provide an estimate of the strength of anharmonicity in a compound. 
One can see from Fig. 7 that the two transverse acoustic branches of phonons have the 
largest Griineisen parameters, indicating that these two branches are most anharmonic and 
interact with the longitudinal acoustic mode strongly. In order to describe quantitatively the 
anharmonic effects, we further calculate the average Griineisen parameter (7) for the three 



acoustic vibration modes by using the similar method proposed in Ref. 38]: 7 = \J {jf). For 
Am0 2 , the values of 7ta, Ita'i Ila are 2.75, 2.48, and 2.16, respectively. And the average 
7 of the three acoustic modes is then calculated to be 2.46. These transverse modes play an 
important role in lattice thermal resistance. Furthermore, since the velocity of transverse 
modes is much lower than the longitudinal one, the lattice thermal conductivity is dominated 
by the lower- velocity transverse modes. 

The acoustic Debye temperature closely related to the lattice thermal conductivity of a 
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FIG. 7: (Color online). Calculated Griineisen parameters for Am02- Blue, red, and olive symbols 
represent for the results corresponding to the two transverse acoustic (TA and TA') and one 
longitudinal acoustic (LA) branches of phonon respectively. The inset depicts the first Brillouin 
zone. 
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FIG. 8: (Color online). Calculated lattice thermal conductivity for A111O2 at different temperatures, 
and corresponding experimental as well as other theoretical results. 



material can be approximately calculated by using the relationship 57 



where G and n are the Debye temperature and number of atoms in the unit cell, respectively. 
For AmC>2, O a is calculated to be 233.0 K. Based on the relationship in Eq. (Q, and the 
above obtained quantities, we calculate the lattice thermal conductivity for Am0 2 and show 
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the results from the temperature of 300 K to 1500 K in Fig. 8. For comparison, previous 



experimental resu 
by Uchida et al. 



ts by Bakker et al. 



59(1 and Nishi et al. |43|, as well as theoretical results 



601 ]. Kurosaki et al. 61[, and Lemehov et al. [62J are also shown in Fig. 



8. Clearly, in the temperature range from 300 to 700 K, our theoretical results accord well 
with the experimental measurements. At the higher temperature range of above 700 K, 
our calculated thermal conductivity are a little smaller than the corresponding experimental 
results. This disagreement might come from two different factors. Firstly, the electronic 
thermal conductivity which occupies larger ratios at higher temperatures is not considered 
here because of limitations of our computational methods. Secondly, the Debye- Callaway 
model we used does not take into account the optical branches of phonons, some of which 
have appreciable group velocities and are able to contribute to thermal conductivity. Over- 
all, the coarse agreement between our theoretical results on thermal conductivity with the 
experimental ones is fortuitous. 



IV. CONCLUSIONS 



In summary, a thorough DFT+[7 study has been performed to investigate the electronic, 
mechanical, and thermodynamical properties of AH1O2. It is found that the chemical bond- 
ing between the americium and oxygen atoms show mixed ionic and covalent characters. The 
valence band maximum and conduction band minimum are contributed by 2p-5f hybridized 
and 5/ electronic states respectively. In comparison with PUO2, the bonding length and 
charge transfer in Am02 are both slightly smaller, indicating stronger covalent and weaker 
ionic interactions in the latter than than in the former. Through calculations of the elastic 
constants and various moduli, we have predicted that compared to PUO2, AH1O2 is more 
stable against shear forces, but less stable against pressures. The stress-strain relationship 
of AH1O2 has also been examined along the three low-index directions by FPCTT calcula- 
tions. Our results show that the [100] and [111] directions are the strongest and weakest 
tensile directions, respectively. In comparison with PUO2, the theoretical tensile strengths 
of Am02 are smaller. At the thermodynamic aspects, the phonon dispersion of Am02 has 
been calculated at different lattice volumes, with the LO-TO splitting effect being modified 
through Born effective charges. Subsequently, we have systematically obtained the heat 
capacities and lattice expansion curve, both of which accord well with the experiments. Fi- 
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nally, by differentiating the phonon dispersion curves at different lattice constants, we have 
further calculated the Gruneisen parameters of Am02 and determined its lattice thermal 
conductance at different temperatures, which show the overall agreement with attainable 
experimental measurements. 
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